The two sided parsec scale structure of the Low Luminosity 
Active Galactic Nucleus in NGC 4278 



M. Giroletti 1 ' 2 , G. B. Taylor 3 ' 4 , and G. Giovannini 1 - 2 



ABSTRACT 

O 

o . 

CN) . We present new Very Long Baseline Interferometry observations of the LINER 

O . galaxy NGC 4278. The observations were taken with the Very Long Baseline Ar- 

ray (VLBA) and a single antenna of the Very Large Array (VLA) at 5 GHz and 

8.4 GHz and have a linear resolution of < 0.1 pc. Our radio data reveal a two 
Os, ~ 

sided structure, with symmetric S"-shaped jets emerging from a flat spectrum 

core. We fit the jet brightness with gaussian components, which we identify from 
a previous observation taken five years before. By comparing the positions of 
! the components in the two epochs, we measure motions between 0.45 ± 0.14 and 

3.76 ± 0.65 mas, corresponding to apparent velocities < 0.2 c, and to ages in the 
Tj- ! range 8.3 — 65.8 years. Assuming that the radio morphology is intrinsically sym- 

metric and its appearance is governed by Doppler beaming effects, we find that 
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NGC 4278 has mildly relativistic jets {(3 ~ 0.75), closely aligned to the line-of- 
sight (2° < 6 < 4°). Alternatively, the source could be oriented at a larger angle 
and asymmetries could be related to the jet interaction with the surrounding 
medium. We also present new simultaneous VLA observations between 1.4 and 
43 GHz, and a 5 GHz light curve between 1972 and 2003. The radio spectrum 
can be fit by a relatively steep power-law (a = 0.54). We find significant variabil- 
ity at 5 GHz. All these arguments indicate that the radiation from NGC 4278 
is emitted via the synchrotron process by relativistic particles accelerated by a 
supermassive black hole. Despite a much lower power, this is the same process 
that takes place in ordinary radio loud AGNs. 

Subject headings: galaxies: active — galaxies: nuclei - - galaxies: individual 
(NGC 4278) 
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1. Introduction 

Although objects hosting an Active Galactic Nucleus (AGN) represent only a small 
fraction of the total number of extragalactic sources, there is growing evidence that a low 
level of nuclear activity is a common feature among galaxies. Objects presenting a spec- 
tral signature of such activity include low-ionization nuclear emission-line region (LINERs, 
Heckman 1980), low luminosity Seyfert galaxies, and "transition nuclei", i.e. nuclei with 
spectra intermediate between LINERs and HII regions. These objects are grouped under 
the name of low luminosity active galactic nuclei (LLAGN, see Ho, Filippenko, & Sargent 
1997). 

Dynamical searches and stellar kinematics have yielded strong evidence for the existence 
of compact dark objects of mass ~ 10 6 to 1O 9,5 M (Kormendy & Richstone 1995; Magorrian 
et al. 1998) in virtually all bulge-dominated galaxies. These are most likely in the form 
of super massive black holes (SMBHs); however, what is still unknown is the state of the 
SMBH, and the reason why they are so often quiescent. In the light of these findings, the 
debate over the physical processes that are responsible for the "activity" in LLAGN has been 
a major issue in recent years. These galaxies may host a scaled down version of the more 
powerful AGN, with an accreting super massive black hole rendered under-luminous because 
of a low accretion rate and/or low radiative efficiency Alternatively, the activity might be 
related to stellar processes (e.g., Shields 1992; Filippenko & Terlevich 1992; Alonso-Herrero 
et al. 2000). 

Along with hard X-ray data (e.g. Ho et al. 2001; Terashima et al. 2002), radio obser- 
vations are a useful tool to discriminate between an AGN and a starburst powered nucleus. 
VLA observations with arcsecond (Filho, Barthel, & Ho 2000, 2002) or subarcsecond (Nagar 
et al. 2000, 2002) resolution have successfully revealed the presence of a compact radio core 
in most LINERs and low luminosity Seyferts, with the non-detections likely caused by limits 
in sensitivity only. The corresponding brightness temperatures are TJ, > 10 2 ' 5-4 K, which are 
consistent with a non-thermal origin. The spectral index of these radio cores is typically flat 
- a well-known signature of synchrotron self-absorbed emission (Nagar, Wilson, & Falcke 
2001). 

High angular resolution Very Long Baseline Interferometry (VLBI) observations yield 
more conclusive evidence: all (16/16) LLAGNs in a flux density and distance limited sample 
possess milliarcsecond compact cores, as revealed with the VLB A at 6 cm (Falcke et al. 2000; 
Nagar et al. 2002). This implies brightness temperatures T B > 10 8 K, which definitely rule 
out any explanation based on thermal processes, and point directly to synchrotron emission 
from the base of a parsec scale jet. Indeed, the five sources with the brightest cores have 
parsec scale jets. Furthermore, high frequency radio spectra can be used to discriminate 
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between pure ADAF (Advection Dominated Accretion Flows, Narayan & Yi 1995) models 
and non-thermal jets (Di Matteo et al. 1999; Di Matteo, Carilli, & Fabian 2001; Anderson, 
Ulvestad, & Ho 2004). 

The nearby elliptical galaxy NGC 4278 has been investigated in detail at most wave- 
bands. In the optical, HST observations reveal a central point source and a large distribution 
of dust located north- northwest of the core (Carollo et al. 1997). Ionized nuclear gas typical 
of a LINER is found in this galaxy (Goudfrooij et al. 1994), possibly associated with an 
external ring of neutral hydrogen in PA 135° (Raimond et al. 1981). 

Radio continuum observations of NGC 4278 on kpc scale have been carried out at 
frequencies between 5 GHz and 43 GHz, revealing a compact source (Di Matteo, Carilli, & 
Fabian 2001; Nagar, Wilson, & Falcke 2001; Nagar et al. 2000). At 1.4 GHz, both the FIRST 
(Becker, White, & Helfand 1995) and the NVSS (Condon et al. 1998) data confirm that the 
source is compact. Only at 8.4 GHz, at a resolution of ~ 200 mas, is NGC 4278 slightly 
resolved into a two-sided source with an extension to the south (Wilkinson et al. 1998). 
Compactness and a flat radio spectrum suggest that the radio emission is non-thermal, and 
in this respect the emission from NGC 4278 seems to be very similar to that of powerful 
radio loud AGNs, such as QSO and BL Lacs; however, the total radio luminosity of the 
source is only Pi.4gh z = 10 216 W Hz -1 , i.e. at least two orders of magnitude less than those 
powerful objects. 

On parsec scales, early VLBI experiments at 18 cm and 6 cm have revealed a core 
dominated structure, with an elongated feature extending to the north-west and possibly to 
the south on scales of some 10 mas (Jones, Terzian, & Sramek 1981; Jones, Sramek, & Terzian 
1982; Schilizzi et al. 1983; Jones, Wrobel, & Shaffer 1984). More recent observations with the 
VLBA at 6 cm have yielded a dramatic improvement in our knowledge of this source: Falcke 
et al. (2000) reveal an extended core and an elongated region to the southeast on scales of 
a few milliarcsecond; Giovannini et al. (2001), thanks to more short spacings provided by a 
VLA antenna in addition to the full VLBA, have detect emission on the opposite side of the 
core as well. However, the discrepancy in the total flux density [(87 ±4) vs. (400 ±20) mJy] 
is more than the simple addition of a VLA antenna can justify; it is possible that the source 
may be variable or simply that there have been problems in the calibration of the flux density 
scale. Bondi et al. (2004) also detect a trace of two-sided emission, although heavily resolved. 
Finally, VLBA phase-referenced observations have succeeded in detecting the source on sub- 
pc scale even at 43 GHz, where it shows only a core and a hint of low level emission to the 
north (Ly, Walker, & Wrobel 2004). 

In the present paper, we consider new VLBA observations at two frequencies (5 GHz 
and 8.4 GHz), taken on 2000 August 27. We compare the new 5 GHz image to previous 



-4- 



images and discuss the morphology and the motion of components. Furthermore, the 8.4 
GHz data provide a favorable combination of resolution and sensitivity, constraining the 
brightness temperature and the location of the central black hole. 

We describe the new observations in §2 and present the results in §3. A discussion is 
given in §4, and the main conclusions are summarized in §5. Throughout this paper, we 
adopt for the Hubble constant a value of 71 km s _1 Mpc" 1 ; however, since NGC 4278 has a 
direct distance measurement of 14.9 Mpc (Tonry et al. 2001; Jensen et al. 2003), we do not 
make use of its redshift (z = 0.00216), which is affected by proper motion in a galaxy group. 
At the distance of NGC 4278, 1 mas corresponds to a linear scale of 0.071 pc. We define the 
spectral index a following the convention that S(u) oc u~ a . 

2. Observations and Data Reduction 

2.1. VLB A observations 

We observed NGC 4278 with an 11 element VLBI array composed of the NRAO Very 
Long Baseline Array (VLBA) and a single 25 m VLA antenna for 10 hours. The observing 
run was performed on 2000 August 27, switching between 5 GHz (with 4 IFs at 4971.49, 
4979.49, 4987.49, 4995.49 MHz) and 8.4 GHz (8405.49, 8413.49, 8421.49, 8429.49 MHz). 
Sixteen channels per IF and full polarization data were recorded. Short scans on 3C279 were 
taken in order to find fringes, and several observations of OQ208 were repeated at different 
hour angles in order to calibrate the polarization leakage terms. J1310+3220 and J1751+0939 
were also observed to perform checks on the amplitude and polarization calibration. 

The correlation was carried out at the AOC in Socorro. The distribution tapes were read 
into the NRAO Astronomical Image Processing System (AIPS) for the initial calibration. 
We followed the same scheme for the data reduction of both 5 GHz and 8.4 GHz data-sets. 
As a first step, we corrected our data using the accurate position information obtained by Ly, 
Walker, & Wrobel (2004) in a phase-referenced experiment with the VLBA at 43 GHz (RA 
12 h 20 m 06 s . 825429, Dec 29° 16' 50".71418). We then performed the usual calibration stages 
(removal of instrumental single band delay, and the R — L delay, and bandpass calibration) 
using scans on 3C279. 

After flagging bad data, we obtained good models for the calibrators, which we used 
to improve the amplitude calibration for the entire data-set. Finally, we removed the in- 
strumental polarization leakage using OQ208. A single-source file for NGC 4278 was then 
created by averaging the frequency channels and IFs and exported into Difmap (Shepherd, 
Pearson, & Taylor 1994, 1995) for editing and self-calibration. 
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Thanks to a good calibration and position information, we could obtain final images 
with only a few iterations of phase self-calibration. One cycle of amplitude self-calibration 
with a long solution interval (30 minutes) has also been performed before obtaining the final 
(u, v)— data. Total intensity images with natural weights are presented in Figs. 1 and 2, and 
the most relevant parameters are summarized in Table 1. Stokes 'Q' and 'IT images were 
also made. 

The final calibrated datasets are not suitable for a direct combination into a spectral 
index map. Different observing frequencies correspond to different coverage of the (u, v)- 
plane and this could introduce spurious results in a spectral index map. Therefore we 
selected an annulus in the (u, v )-plane, consisting of only baselines in the range 0.9 MA — 145 
MA. We also tapered the data at 100 MA in u and 60 MA in v in order to account for 
the more elliptical coverage of the 5 GHz data. After this processing, the images have a 
resolution of 3.2 mas x 1.9 mas in PA —2° and were combined to produce the spectral index 
map using the AIPS task COMB. 

We also re-analyzed VLBA+Y1 5 GHz data obtained in 1995 (July 22), taking advantage 
of the new position (Ly, Walker, & Wrobel 2004). The position used in the original correlation 
was off by more than one arcsecond from the actual position of the source. This resulted in 
exceedingly large delays which caused a lot of failures in the fringe fitting process. The new 
accurate positional information allows us to recover a large fraction of useful data, resulting 
in a significantly improved, high fidelity image. We present and briefly discuss the results 
from the reprocessed data in Sect. 3.2. 

2.2. VLA observations 

New VLA observations were taken at 1.4, 5.0, 8.4, 15, 22, and 43 GHz in A configuration 
on 2003 August 17 and reduced in AIPS following standard procedures. 3C286 was used to 
set the flux density scale. Phase referencing to the nearby calibrator 1221+282 failed at high 
frequencies due to bad weather (rain). Thus, we only have an upper limit at 43 GHz and 
the value at 22 GHz is more uncertain than at lower frequencies - we estimate a fractional 
error of 2% up to 8.4 GHz, of 5% at 15 GHz, and of 10% at 22 GHz. 
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3. Results 

3.1. VLBA on 27 August 2000 

The final images reveal a source dominated by a central compact component, with 
emission coming from either side. To the southeast, a jet-like feature extends for ~ 6.5 
mas in PA 155° (measured north to east), then progressively bends into PA 100°. In total, 
the jet is almost 20 mas long, which corresponds to ~ 1.4 pc. On the opposite side, the 
main component is slightly elongated to the north in the 5 GHz map, and the 8.4 GHz data 
clearly show a secondary component in PA —40°. Then, this jet-like feature bends to the 
west turning into a diffuse, uncollimated, low brightness region. 

In total, the source extends over ~ 45 milliarcsecond, i.e. about 3 parsecs. Given 
its morphology and dimension, NGC 4278 meets all requirements for the classification as 
a Compact Symmetric Object (CSO), except that it is under-luminous compared to most 
known CSOs (Readhead et al. 1994; Peck & Taylor 2000). 

The total flux density measured in our images is 120 mJy at 5 GHz and 95 mJy at 
8.4 GHz. If we compare these values with those obtained with the VLA (162 mJy and 114 
mJy, respectively, see Table 5), we find about a 20-25% difference, which can be ascribed 
to the VLBA resolving out some extended emission, probably in the western region. The 
monochromatic luminosity at 1.4 GHz is 3.18 x 10 21 W Hz" 1 , and 2.52 x 10 21 W Hz" 1 at 8.4 
GHz. 

The visibility data are well fit by a five component model at both frequencies. The 
position and dimension of the components are illustrated in Fig. 3; our choice for labeling 
the components is based on their most likely epoch of ejection, as discussed in § 3.3. Table 2 
shows the relevant parameters for all components, re-referenced to the central component 
C. To compute positional uncertainties for the components, we repeated the modelfit on the 
'LL' and 'RR' correlations independently; the semi-dispersion of the results is taken as a la 
error on the position of each component. This yields Ar = 0.04 mas for components .S2 and 
SI, 0.01 mas for N3 and 0.34 mas for N2. In Columns 7 and 8 we report the flux density at 
5 and 8.4 GHz respectively, and in Column 9 we compute the spectral index. Besides being 
the most compact feature, component C presents also the flattest spectral index (a = 0.2) 
and its identification with the core is straightforward. 

Finally, no polarized signal is detected from our data. Both the images in Stokes 'Q' and 
'IT are purely noise-like, with on source la noise levels of 66 and 63 /iJy/beam at 5 GHz and 
of 62 and 60 mJy/beam at 8.4 GHz, respectively. Based on the on-source noise in the 'Q' 
and 'IT images, we place a limit on any polarized signal for the source at < 125 /iJy/beam 
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at 5 GHz and < 120 /zJy/beam at 8.4 GHz (3a values). At the position of the brightest 
component in the total intensity images, this corresponds to a fractional polarization of 
< 0.3% at 5 GHz and < 0.5% at 8.4 GHz. 



3.2. VLBA on 22 July 1995 

The reprocessing of the July 1995 data has yielded a significant improvement over the 
previously published results. The new image (Fig. 4) is in good agreement with the data 
taken at 5 GHz in 2000, showing a central compact component and jet-like features in 
the same position angle. The overall morphology of the source is similar to that shown 
in Giovannini et al. (2001) and we confirm the reality of the emission coming from the 
northwest. However, the noise level has improved by over a factor of 6 (65 vs 400 /iJy) and 
the flux density scale is more consistent with both our new observations as well as the results 
presented in Falcke et al. (2000) and Bondi et al. (2004). The resulting total flux density is 
135 mJy. 

The same five component model that fits the 2000 epoch data has been applied to this 
data set as well, allowing for the components to change in flux density and position. The 
resulting fit is good (reduced x-squared of 0.98) and is presented in Table 3. Uncertainties 
are computed in the same way as for epoch 2000.652, and are 0.13, 0.11, 0.09, and 0.56 mas 
for S2 SI, A3, and A2, respectively. 

An independent fit with only four components has also been tried, and fits the data 
almost as well as the five component model (reduced x-squared of 1.00): components S2, 
SI, and N2 are maintained in both models, while C and A3 can actually be fitted with one 
Gaussian only (C). However, this model gives a less accurate representation of the core 
region, with the component C too extended (3.2 mas x 1.3 mas) to match the size of the 
compact core. We also note that if A3 is removed from the model for the 2000.652 epoch 
data, any following iteration of the fitting algorithm will reintroduce it at the expense of 5*2. 

3.3. Component motion and variability 

If we compare data taken at the same frequency in different epochs, we can get informa- 
tion on the evolution of the source. A first comparison can be made by directly overlaying 
the contours at the same frequency from the two epochs (Fig. 5). Before overlaying them, 
we have convolved the two images with the same restoring beam. We also plot the same 
contours, regardless of the lower noise level of the 2000.652 image. This guarantees that 
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any difference is due to a real change in the source structure and not simply to the better 
sensitivity of one observation. 

The overlay is suggestive of a displacement of the northwest low level emission, mov- 
ing away from the core. Minor differences can be seen also in the core region and in the 
southeastern jet. In general, the blue 2000.652 contours trace the emission in a region more 
extended than the red 1995.551 ones. This behavior can be better quantified if we consider 
the relative motion of the components identified by the modelfit process. Taking the posi- 
tion of the core (component C) as a reference, and assuming it is fixed, we have compared 
the position of the other components. We report the results in Table 4: column (1) labels 
the components, columns (2) and (3) give the relative motion over the five years in polar 
coordinates (distance in mas, angle in degrees); columns (4) and (5) report the apparent 
velocity in mas yr _1 and units of c, respectively; finally, the corresponding kinematic age is 
given in column (6). The radial distance of each component has increased over the five years 
lag between the observations. The motion is larger on the northwestern side; in particular, 
the largest displacement is found for component N2. The corresponding apparent velocities 
are given in Column (5), with the fastest component moving at ~ 0.17 c. 

Assuming that the apparent velocity is constant for each component, we derive ages as 
reported in column (6). S2 and N2 have ages that are consistent, and they could have been 
ejected together about 25 years before epoch 2000.652. SI is the oldest component, and its 
counterpart in the main jet is not detected, probably because it is too distant and extended. 
Finally, N3 is the youngest component, ejected only three years before our first epoch of 
observations; it is likely that a corresponding component S3 has emerged in the counterjet 
but is still confused with the core. Note that the core is the only component whose flux 
density is larger in 2000.652 than in 1995.551. 

Finally, we can compare the total flux density detected in the two epochs. Our data 
indicate a variation of the total flux density of about 10%, with the source being less powerful 
in the second epoch. Note that Falcke et al. (2000) report a value of 87.3 mJy; however, 
that observation was performed without the addition of a VLA antenna, which results in a 
lack of short spacings sensitive to low level emission. As a consequence, component N2 and 
(partly) SI are resolved out, accounting for the difference in total flux density. 

3.4. Spectral index 

We now compare data taken at the same epoch but with different frequencies, in par- 
ticular we present the map of the spectral index from 2000 August 27 (Fig. 6). The features 
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of the spectral index map are in agreement with the single component values determined by 
the model-fit. The core region presents the flattest spectrum (a ~ 0.0), steepening along the 
southern jet (a > 0.5). An irregular spectral index pattern is found in regions with lower 
level emission, such as the outermost edge of the southern jet and the extended northern 
component N2. However, we can still determine an average spectrum for this region, which 
has a = 0.8. 

In Fig. 7 we also plot the simultaneous spectrum taken with the VLA in A configuration 
from 1.4 GHz to 43 GHz. The measured flux densities are reported in Table 5. The best-fit 
power-law spectral index is a = 0.54. 

4. Discussion 

4.1. Jet orientation and velocity 

Our images detect low level emission to the northwest of NGC 4278 with unprecedented 
resolution and sensitivity. We find, in the inner part of the jet, a compact component N3, 
and at a larger distance, ~ 10 mJy of flux density in the region N2. Thus, we classify 
NGC 4278 as a two-sided source, similarly to a few other LLAGN previously studied, e.g. 
NGC 4552 (Nagar et al. 2002), NGC 6500 (Falcke et al. 2000), and NGC 3894 (Taylor, 
Wrobel, & Vermeulen 1998). 

We are therefore faced with the puzzle of determining which of the jets is pointing toward 
us and the direction of the source main axis with respect to the line-of-sight. Although the 
southern jet looks more collimated, the total flux density is larger in the northern components 
than in the southern ones. The high resolution 8.4 GHz image (Fig. 2) clearly shows that 
the inner jet is brighter to the north than to the south. The apparent motion of the northern 
components are also larger than those of the southern ones. Finally, the 43 GHz image in 
Ly, Walker, & Wrobel (2004) shows a hint of emission only to the northwest of the core on 
scales < 1 mas. We therefore argue that the main, approaching jet is the northern one. 

4-1.1. Jet Asymmetries from Doppler Boosting 

To estimate the orientation 9 and intrinsic velocity (3 of the jet, we initially consider a 
simple beaming model, which assumes that pairs of components are ejected simultaneously 
from the core with the same intrinsic velocity and brightness. We apply this model to the 
component pair N2/S2, which has been ejected simultaneously according to our motion 
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measurements (§ 3.3), and we further assume that N2 is the one on the approaching side. 
In this model, the ratio between the arm length r and the proper motion \i of the two 
components are related by 

^S2 r S 2 1 - P cos 9 

From our modelfits we derive that r N2 /rs2 = 7.2 ± 0.2 and HN2I ^S2 = 8 ± 3, so we 
estimate that 4 < R < 10. The arm length ratio corresponds to selecting the hatched area 
between the two solid lines in the (/3, 6?)-plane shown in Fig. 8. The dot-dash lines represent 
the possible combination of (3 and 9 resulting from the apparent separation velocity of the 
two components, which is expressed by the relation (5 sep = (2/3 sin 9)/(l — [3 2 cos 2 9). Finally, 
since we measure motion on both sides and we know the source distance, we can directly 
solve for 9 and f3 (see e.g. Mirabel & Rodriguez 1994); this corresponds to the dashed ellipse 
centered on 9 = 2.7°, (3 = 0.79. 

In principle, one could also consider the brightness ratio between the two components; 
however, as an effect of relativistic time dilation, we are watching the components at different 
stages of evolution. Since we know little about the time evolution of the jet components, 
this hinders the application of the brightness ratio constraint; in any case, a check that 
Sn2/Ss2 > 1 is possible and is consistent with our interpretation. 

Based on the above analysis (Fig. 8), we find mildly relativistic velocities of f3 ~ 0.76 
(r ~ 1.5), and an orientation close to the line-of-sight (2° < 9 < 4°). The resultant Doppler 
factor is 5 ~ 2.7 and the small viewing angle explains the bends visible in both jets as the 
amplification caused by projection effects of intrinsically small deviations, which are common 
in high power radio sources. 

4-1.2. Jet Asymmetries from Interactions 

As an alternate explanation to that proposed above, it is possible that the true radio 
morphology of NGC 4278 is influenced by more than relativistic beaming. For example, there 
may be significant interactions between the jet and the surrounding medium. This could 
explain the lower velocity for the components in the southern side of the jet and cause the 
flux density of SI to be enhanced. If this is the case, the task of identifying corresponding 
pairs of components either side of the core becomes confused. The radio source could then 
be oriented at a larger angle to the line-of-sight than we deduce from beaming arguments. 
In this scenario, the large bends visible in our images would be real and not amplified by 
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geometrical effects, suggesting a strong interaction of the jets with the surrounding medium. 
The combination of a strong interaction with a dense medium and a low power core could 
well account for the small size of NGC 4278. It is unlikely, in any case, that more than 
mildly relativistic velocities are present in the jet of NGC 4278. 

It is interesting to note that Nagar et al. (2002) have detected compact cores in LLAGNs 
almost exclusively in Type 1 objects. Among the possible explanations, they invoke the 
unified scheme, assuming that all LLAGNs have accreting black holes and that we only 
see radio emission in Type 1 objects because of beaming and/or obscuration. Our resulting 
viewing angle determined in §4.1.1 seems to support this picture; however, the Lorentz factor 
is not as high as in more powerful classic radio galaxies, where 3 < T < 10 (Giovannini et 
al. 2001). Therefore, a source similar to NGC 4278 but oriented at 60°, would have a similar 
Doppler factor S ~ 1.1, hence still be detectable. We note that if LLAGNs have low Lorentz 
factor (i.e. T < 3), the Doppler factor does not depend strongly on the viewing angle and 
we should be able to see also sources oriented at large 8. The lack of raido cores in Type 2 
LLAGN then may require alternative explanations such as free-free absorption by a torus. 



4.2. History of emission 

From the result of the modelfit, small to moderate apparent velocities (< 0.2 c) are 
found for the four jet components. Under the assumption of constant velocity, we derive 
that these components must have been ejected from the core between 8.3 and 65.8 years 
before epoch 2000.652 (see Column [6] in Table 4). 

However, these jet components are not to be confused with the hot spots demarcating 
the end of the jet as found in more powerful CSOs (Owsianik & Conway 1998; Peck & Taylor 
2000; Giroletti et al. 2003). Therefore, a kinematic estimate of the real age of the source 
is not possible, and that of SI can only be taken as a lower limit. The low brightness and 
large size of N2, as well as the non-detection of Nl, suggest that components are continually 
ejected from the core, but soon disrupt, without being able to travel long distances and form 
kiloparsec-scale lobes. 

We do not expect that this source will evolve into a kiloparsec-scale radio galaxy, but 
rather that it will only periodically inflate slowly, as visible in Fig. 5. The relatively low 
velocity jets discussed in § 4.1 cannot bore through the local ISM and escape, as demonstrated 
by the lack of hot spots which are prevalent in higher power CSOs. This behavior can be 
ascribed to a low power central engine, which cannot create highly relativistic jets. 

In Fig. 9 we plot the flux density history for NGC 4278 at 6 cm using data taken at 
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the Westerbork Synthesis Radio Telescope, the Green Bank 300 ft radio telescope and the 
VLA. A previous light curve was published by Wrobel & Heeschen (1991), to which we add 
12 points, from observations obtained between 1972 and 2003. The light curve shows that 
the radio source is variable, prone to both outbursts and low states. A burst is certainly 
present around 1985, while in more recent years the source has been showing less activity. 

It is difficult to connect the burst with the ejection of new components, both because of 
the uncertainties related to the age of a single component, and the possible time lag between 
component ejection and total flux density enhancement. It is clear however that the source 
presents a high degree of variability, most likely related to the presence of an active nucleus. 
This is an interesting result, suggesting that repeated observations of other LLAGN are 
desirable, given that few studies on this topic exist so far (Wrobel 2000; Nagar et al. 2002). 

4.3. The origin of energy in LLAGN 

Thanks to its proximity, NGC 4278 makes an excellent laboratory to study the properties 
of LLAGNs and to explore the nature of the physical processes taking place. In particular, 
there are three important clues for solving the puzzle over the origin of the radio emission: 
(1) the shape of the radio spectrum; (2) the parsec-scale morphology, and (3) the brightness 
temperature of the VLBI core. 

In 7 nearby galaxies studied by Di Matteo et al. (1999); Di Matteo, Carilli, & Fabian 
(2001), ADAF models are unsuccessful in fitting the high frequency radio spectra, and syn- 
chrotron radiation from non-thermal particles in jets must contribute significantly. Con- 
versely, Anderson, Ulvestad, & Ho (2004) have revealed in six LLAGN VLBA spectral in- 
dexes that are remarkably flat or even slowly rising up to 43 GHz. Although these cores 
are still too luminous to be explained in terms of only radiation from an ADAF, the jet 
properties remain to be determined and successfully imaged in these sources. 

In our case, the overall VLA spectrum can be fit by a power-law with an intermediate 
index a = 0.54 (see §3.4). A possible steepening occurs above 22 GHz, which argues against 
the presence of an ADAF as significantly contributing to the total emissivity. A low frequency 
flattening/turnover may be present around 408 MHz, with S408 mh z — 650 mJy (Colla et al. 
1972) and S365 mh z — 569 mJy (Douglas et al. 1996). Variability and resolution effects are 
also important in the study of the spectrum of this source and other AGN - for example, the 
flux density from the low brightness northwestern region was missed by Falcke et al. (2000), 
resulting in a flatter high-frequency spectral index. In any case, the radio spectrum between 
151 MHz and 43 GHz is globally similar to that of more powerful radio loud AGN, i.e. to 
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that of synchrotron emission from relativistic particles in a jet. 

At higher frequencies, NGC 4278 possesses an optical central compact core which lies 
nicely at the low power end of the correlation between optical and radio core luminosity in 
FR I radio galaxies (Capetti et al. 2002). This correlation can be interpreted in terms of a 
common non-thermal synchrotron origin, with the radiation being produced in relativistic 
jets (Capetti & Celotti 1999). 

The parsec-scale morphology revealed by our VLBA images is clearly in agreement 
with this scenario, showing twin jets. This result and the high brightness temperature 
Tb = 1.5 x 10 9 K are in agreement with similar findings from Nagar et al. (2002); therefore, 
we conclude that the emission from NGC 4278 and from similar LLAGNs originates in radio 
jets via the synchrotron process. 

In this scenario, detection of polarized signals would be an interesting signature of the 
synchrotron process; indeed, VLA observations at 8.4 GHz have found weak linear polar- 
ization in NGC 4278 at a level of 0.34% and in other LLAGN, although the detection rates 
and percentage of linear polarization in these sources are lower than in more powerful AGNs 
(Bower, Falcke, & Mellon 2002). Our new polarization data for NGC 4278 show that it is 
unpolarized on milliarcsecond scales, as also found by Bondi et al. (2004). This could be 
due to limits in sensitivity. Since synchrotron radiation is expected to have an intrinsically 
high percentage of linear polarization, one needs to invoke a significant amount of Fara- 
day depolarization, or some other effect, to explain the low amount of linear polarization 
observed. Note that low or absent linear polarization is a characteristic of CSOs (Peck & 
Taylor 2000; Giroletti et al. 2003), although recent observations have found exceptions to 
this rule (Gugliucci et al. 2004). 

Finally, it is interesting to discuss the position of the central black hole responsible for 
the accretion and launch of the jets. Our observations allow us to image the nuclear region 
of NGC 4278 on scales of < 0.1 pc (1 mas = 0.071 pc). Based on the correlation between 
the black hole mass and the central velocity dispersion, we can estimate Mbh for NGC 4278. 
The relation takes the form log(MBH/-M ) = 8.13 + 4.02 log(cr/cr ), where o"o = 200 km s _1 
(Tremaine et al. 2002) and o = (258 ± 15) km s" 1 (Barth, Ho, & Sargent 2002). It follows 
that M B h = 3.7 x 1O 8 M and that the Schwarzschild radius 1R S = 3.5 x 10~ 5 pc; therefore, 
our images probe scales of some 10 3 i?s. 
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5. Conclusions 

We present new VLBA and VLA data for the nearby (d = 14.9 Mpc) LLAGN NGC 4278. 
Our new VLBA data show a two-sided emission on sub-parsec scales in the form of twin jets 
emerging from a central compact component (Tg = 1.5 x 10 9 K), in a similar way to that seen 
in more powerful radio loud AGNs. In agreement with this, the spectral index distribution 
for the radio source reveals a flat spectrum core region with steep spectrum jets on either 
side. 

By comparison with previous observations, we discover proper motion for components 
in both jets over a five years time baseline. We find low apparent velocities (< 0.2 c) for 
the jet components and estimate the epoch of their ejection as 10 — 100 years prior to our 
observations. Based on our analysis, we suggest that the north-west side is the approaching 
side, and that the jets of NGC 4278 are mildly relativistic with (3 ~ 0.75. 

The central black hole in NGC 4278 is active and able to produce jets, which are 
responsible for the bulk of the emission at radio to optical frequencies in this LLAGN. 
However, the lifetime of the components of < 100 years at the present epoch, combined 
with the lack of large scale emission, suggests that the jets are disrupted before they reach 
kiloparsec scales. 

The study of the flux density history at 6 cm between 1972 and 2003 shows significant 
variability (> 100%) on time scales of a few years, which might be related to the ejection 
of new components. This subject needs to be explored for other LLAGNs as well, as it can 
give insights into the state of the central black hole in these sources. 
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Table 1. Image Parameters 







beam 


beam 






peak 


total 




Epoch 


Freq. 


size 


P.A. 


noise 




brightness 


flux density 


P 




(GHz) 


(mas x mas) 


(°) 


(/iJy beam" 




(mJy beam -1 ) 


(mJy) 


(W Hz -1 ) 


1995.551 


5.0 


3.0 x 1.8 


-25 




65 


49.5 


135 ±7 


3.57 x 10 21 


2000.652 


5.0 


3.0 x 1.7 


-7 




45 


40.5 


120 ±6 


3.18 x 10 21 


2000.652 


8.4 


1.8 x 1.0 


-4 




50 


24.6 


95 ±5 


2.52 x 10 21 



Table 2. Results of the modelfit to epoch 2000.652 



Component 


r 


e 


a 


b/a 


$ 


5*5 GHz 


58.4 GHz 


a 




(mas) 


(°) 


(mas) 




(°) 


(mJy) 


(mJy) 




C 


0.00 


0.0 


0.99 


0.48 


-26.6 


31.5 


28.0 


0.2 


S2 


2.57 


162.0 


1.65 


0.36 


76.4 


7.5 


3.0 


1.8 


SI 


8.51 


142.8 


7.39 


0.50 


-68.3 


28.0 


20.8 


0.6 


N3 


1.96 


-25.1 


3.23 


0.74 


-39.0 


45.4 


36.6 


0.4 


N2 


18.46 


-67.7 


11.17 


0.45 


-71.9 


8.7 


6.3 


0.6 



Note. - Parameters of each Gaussian component of the model brightness distribution: 
r, 9, polar coordinates of the center of the component relative to an arbitrary origin, with 
polar angle measured from north through east; a, b, major and minor axes of the FWHM con- 
tour; <3>, position angle of the major axis measured from north through east; S5GHZ, S8.4GHz, 
flux densities; a, spectral index between 5 and 8.4 GHz. 
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Table 3. Results of the modelfit to epoch 1995.551 



Component 


r 


e 


a 


b/a 


$ 


5*5 GHz 




(mas) 


(°) 


(mas) 




(°) 


(mJy) 


C 


0.00 


0.0 


0.99 


0.48 


-26.6 


26.7 


S2 


2.13 


164.8 


1.65 


0.36 


76.4 


10.8 


SI 


7.85 


142.6 


6.00 


0.58 


-71.1 


32.8 


N3 


1.32 


12.1 


2.86 


0.46 


-23.4 


52.5 


N2 


14.80 


-64.7 


12.51 


0.32 


-67.5 


13.5 



Note. - Parameters of each Gaussian component of the model brightness distribution: 
r, 6, polar coordinates of the center of the component relative to an arbitrary origin, with 
polar angle measured from north through east; a, b, major and minor axes of the FWHM 
contour; $, position angle of the major axis measured from north through east; Ssghz, flux 
density at 5 GHz. 



Table 4. Component motion at 5 GHz 



Component 


Ar 


P 


fi 


/^app 


age 




(mas) 


(°) 


(mas yr _1 ) 


Kpp/C) 


(yrs) 


C 


reference 










S2 


0.45 ±0.14 


148.7 


0.088 ±0.027 


0.020 ±0.006 


29.1 ±9.3 


SI 


0.66 ±0.12 


145.1 


0.129 ±0.024 


0.030 ±0.006 


65.8 ± 12.4 


N3 


1.21 ±0.09 


-66.3 


0.237 ±0.018 


0.055 ±0.004 


8.3 ±0.5 


N2 


3.76 ±0.65 


-79.5 


0.737 ±0.128 


0.171 ±0.030 


25.0 ±4.8 
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Table 5. VLA Flux Densities 



Frequency 


Flux Density 


(GHz) 


(mJy) 


1.4 


331 ±7 


5.0 


162 ±3 


8.4 


114 ±2 


15 


99 ±5 


22 


73 ±7 


43 


< 50 
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NGC 4278 5 GHz August 2000 
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Fig. 1.- 
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is from — 



VLBA+Y1 image of NGC 4278 at 5 GHz, 2000 August 27. Contours are drawn at 
.., 128) times the lowest contour, which is 0.150 mJy/beam. The grey scale range 
0.2 to 40.5 mJy/beam. 
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Fig. 2.— VLBA+Y1 image of NGC 4278 at 8.4 GHz, 2000 August 27. Contours are drawn 
at (1, 2, 4, 128) times the lowest contour, which is 0.150 mJy/beam. The grey scale range 
is from —0.2 to 24.6 mJy/beam. 
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Fig. 3. — Model components for epoch 2000.652, overlaid on the lowest contour from the 5 
GHz image (150 /xJy/beam). Parameters of the components are given in Table 2. 
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NGC4278 5 GHz July 1995 
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Fig. 4- 
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VLBA+Y1 image of NGC 4278 at 5 GHz, 1995 July 22. Contours are drawn at 
.., 128) times the lowest contour, which is 0.195 mJy/beam. The grey scale range 
0.3 to 49.5 mJy/beam. 
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Fig. 5. — An overlay of the intensity contours at 5 GHz for epoch 1995.551 (red) and 2000.652 
(blue). Both images have been convolved with the 1995.551 beam (3.0 mas x 1.8 mas in PA 
—25°) and plotted with the same contours, i.e. (1, 2, 4, ...) x 0.2 mJy/beam 
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Fig. 6. — The spectral index map between 5 GHz and 8.4 GHz with overlaid contours at 
5 GHz. The color scale range is between —0.8 (red, inverted indices) and 1.5 (blue, steep 
spectrum); contours are drawn at (1, 3, 10, 30, 100) times the lowest contour, which is 0.2 
mJy/beam. In our convention, S(v) oc v~ a . 
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Fig. 7. — VLA flux density measurements between 1.4 GHz and 43 GHz. 
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Viewing angle 8 [°] 

Fig. 8.— (0,/3)-plane for NGC 4278, as derived from the component pair N2/S2. The 
solid lines are the limits from the arm length ratio (r N2 /rs2 = 7.2 ± 0.2, see discussion), 
while the dot-dash lines represent the results from the apparent separation velocity ((3 scp = 
0.191 ± 0.030); finally, the motion ratio and the measured linear distance of the source select 
the region delimited by the dashed ellipse. 
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Fig. 9. — The light curve for NGC 4278 at 6 cm. The asterisks mark the estimated epoch of 
emission of components S2, N2, and iV3. The SI ejection is too far back in time to appear 
in this plot. Different symbols mark data taken with the WSRT (filled triangles; see Ekers, 
Fanti, & Miley 1983; Schilizzi et al. 1983), the VLA (filled squares; Wrobel & Heeschen 
(1991) and references therein; Nagar, Wilson, & Falcke (2001); present work), and the Green 
Bank 300 ft telescope (empty pentagons; Sramek 1975; Langston et al. 1990; Becker, White, 
& Edwards 1991; Gregory & Condon 1991; Gregory et al. 1996). 



